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Abstract 

A method for extracting time-varying oscillatory mo- 
tions from time series records is applied to Lagrangian 
trajectories from a numerical model of eddies generated 
by an unstable equivalent barotropic jet on a beta plane. 
An oscillation in a Lagrangian trajectory is represented 
mathematically as the signal traced out as a particle or- 
bits a time-varying ellipse, a model which captures wave- 
like motions as well as the displacement signal of a par- 
ticle trapped in an evolving vortex. Such oscillatory fea- 
tures can be separated from the turbulent background 
flow through an analysis founded upon a complex- valued 
wavelet transform of the trajectory. Application of the 
method to a set of one hundred modeled trajectories 
shows that the oscillatory motions of Lagrangian par- 
ticles orbiting vortex cores appear to be extracted very 
well by the method, which depends upon only a hand- 
ful of free parameters and which requires no operator 
intervention. Furthermore, vortex motions are clearly 
distinguished from wavelike meandering of the jet — the 
former are high frequency, nearly circular signals, while 
the latter are linear in polarization and at much lower 
frequencies. This suggests that the proposed method can 
be useful for identifying and studying vortex and wave 
properties in large Lagrangian datasets. In particular, 
the eccentricity oi the oscillatory displacement signals, a 
quantity which is not normally considered in Lagrangian 
studies, emerges as an informative diagnostic for charac- 
terizing qualitatively different types of motion. 

1 Introduction 

Understanding the role of long-lived eddies in the global 
ocean circulation is a major topic in oceanographic 
research. Lagrangian floats and drifters constitute 
an invaluable platform for observing the relatively 
small spatial and temporal scales associated with such 
vortices, and indeed many dozens of publications have 
examined vortex properties from regional Lagrangian 
experiments, e.g. Armi et al. (1989), {Flament et ai, 
2001), and {Shoosmith et ai, 2005). Systematic investi- 
gation of the now-extensive historical set of Lagrangian 
drifter and float trajectories is, however, hampered 
by a technical limitation — the lack of a reliable and 
precise method to meaningfully separate vortex motions 
from the background flow. Various approaches have 
been proposed (e.g. Armi et at, 1989; Flament et al, 
2001; Lankhorst, 2006). However, the fact that this 
problem remains unsolved is evidenced by the fact that 
large-scale studies either continue to rely on traditional 
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subjective identification (e.g. Shoosmith et ai, 2005), 
or else to focus on measures of the effect of vortices on 
trajectories rather than the properties of the vortices 
themselves [Griff a et al, 2008). 

Recently a new and general method has been de- 
veloped [Lilly and Gascard, 2006; Lilly and Olhede, 
2009a, 2011), grounded in nonstationary time series 
theory, which permits the automated identification, 
extraction, and analysis of time-varying oscillatory 
features of unknown frequency — such as the signature 
of a particle trapped in a vortex or advected by a 
wave. Here the method is applied to Lagrangian 
trajectories from a numerical simulation in order to 
illustrate the possibility of accurately extracting and 
distinguishing vortex currents and wavelike motions. 
All relevant analysis software is freely distributed to the 
community as a part of a Matlab(R) toolbox, available 
at http://www.jmlilly.net. 

2 Numerical Model 

For an idealized numerical model generating eddies as 
well as background variability, we choose an equivalent 
barotropic quasigeostrophic model of an initially un- 
stable jet on a beta plane. The model integrates the 
equation for conservation of potential vorticity following 
the geostrophic fiow 
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where $ is the streamfunction, /3 = df/d-dis the 
derivative of the Coriolis frequency / with latitude 
the parameter L^i is the Rossby radius of deformation, 
and k is the vertical unit vector. The model is initialized 
at time t = Q with an eastward jet of strength U and 
width Y having a profile, with i being the eastward unit 
vector, given by 
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which corresponds to a maximum initial vorticity 
anomaly within the jet of Co = ttU /[2Y). 

Parameters are chosen to give a strong jet with a defor- 
mation radius that is small compared to the radius of the 
earth. The central latitude [y = 0) at the jet axis is set 
to d = 45° N, the jet width 2Y and deformation radius 
Lo are both 80 km, and the maximum initial velocity 
is 2.08 m s~'^ . These choices give a jet Rossby number 
Ro = C,o/f = 0.80, and a value of /3Lr,/Co of l/(207r) » 
1/60 — meaning that the jet relative vorticity anomaly is 
much larger than the change in planetary vorticity over 
one deformation radius, or over the jet width. This sys- 
tem is a convenient way of generating eddies, and is not 
intended to represent a particular oceanic current. The 



1 



initial condition is let to freely evolve for 360 days with a 
time step of 5 x 10^^ x 360 days « 26 minutes in a domain 
of length 27r x 400 « 2500 km on each side. The model 
is seeded with 100 drifters that all initially lie along the 
y = line, sampled every 10 time steps or 4.3 hours. 

A snapshot of the model is shown in Fig. 1, together 
with overlays of ellipses characterizing oscillatory 
Lagrangian variability from our subsequent analysis. 
Vortices are generated by barotropic instability and 
then tend to drift westward as well as meridionally 
due to nonlinear beta drift (e.g. Lam and Dritschel, 
2001), with anticyclones propagating equatorward and 
cyclones propagating poleward. Dipole interactions 
are also seen, as captured by the red/blue pair of 
ellipses in the upper left quadrant of Fig. 1, although 
these tend to eventually break down. An informative 
animation of snapshots such as the one shown in 
Fig. 1 over the entire model run duration, but with the 
ellipses color coded by drifter number for visual clarity, 
may be found in the file vortexmovie .mov at ftp: 
//f tp .nwra. com/outgo ing/lilly/vortex. Trajecto- 
ries are shown in Fig. 2a, with the dominating presence 
of vortex motions apparent as tightly looping or cycloidal 
curves reaching northwestward and southwestward. 

3 Analysis Method 

A modulated elliptical signal, introduced by Lilly and 
Olhede (2010), is a time- varying oscillation in two 
dimensions. Such a signal is expressed in matrix form as 



Snapshot of Unstable Jet with Lagrangian Ellipses 



cos9{t) ~sme{t) 
sm0{t) cos 9{t) 



a{t) cos <j){t) 
b{t) sin0(i) 



(3) 



which is the parametric equation for a particle orbiting 
a time-varying ellipse with semi-major and semi-minor 
axes a{t) and \b{t)\, where a{t) > \b{t)\ > 0, and with its 
major axis inclined at an angle 6{t) with respect to the 
X-axis. The phase (/)(t) gives the instantaneous position 
of the particle along the periphery of ellipse. The 
particle orbits the ellipse in the mathematically positive 
or negative direction according to the sign of b{t). 

Details of the modulated elliptical signal, including 
conditions for associating a unique set of time-varying 
ellipse parameters to a given oscillatory signal x(i), 
are discussed in Lilly and Olhede (2010). An impor- 
tant special case is that of a familiar pure sinusoidal 
oscillation in two dimensions; however, the model (3) 
is considerably more general. A practical constraint is 
that the ellipse properties should vary slowly compared 
with the timescale 2t: / ■^(j){t) over which the particle 
orbits the ellipse, in order that the subsequent analysis 
method have small errors [Lilly and Olhede, 2011). 

The modulated elliptical signal with slowly-varying 
ellipse parameters is a good model for the displacement 
signal of a particle trapped in an evolving vortex. This 
class of signals includes Lagrangian displacements due 
to steady circular vortex solutions, steadily strained or 
sheared elliptical anticyclones (Ruddick, 1987), and low- 
frequency periodic oscillations of an elliptical shallow wa- 
ter vortex {Young, 1986; Holm, 1991), all observed with 
instruments that may be experiencing a drift through 
the vortex in addition to the vortex currents themselves. 

Ellipse size, shape, and frequency are usefully char- 
acterized as follows. The ellipse shape is described by 




Figure 1: A snapshot of an unstable eastward equivalent 
barotropic jet on a mid-latitude beta plane at day 
240, with domain details as discussed in the text. The 
shading is the absolute value of the quasigeostrophic 
potential vorticity appearing in (1). Estimated instan- 
taneous ellipses due to Lagrangian oscillatory motions, 
created as in Section 3, are overlaid. Highly eccentric 
ellipses with eccentricity e > 0.95 are shown in green, 
while positively-rotating and negatively-rotating ellipses 
with e < 0.95 are shown in red and blue, respectively. 
Dots mark the instantaneous locations of 100 La- 
grangian particles initially deployed along the jet axis. 
Black dots indicate particles in which an oscillatory 
signal is detected at this moment, while white dots 
indicate the locations of other particles. 



the eccentricity e{t) = \J\ — b'^{t)j a^(t). The geometric 
mean radius and geometric mean velocity are defined as 



R(t) = ^/ a(t)\b{t)\, V(i) 



sgn{6u(t)} Vau(i)|6u(t)| 
(4) 

where the latter quantity is found by writing the velocity 
u(t) = 2jx(i) as a time-varying ellipse of the form (3) 
but with a different set of ellipse parameters, denoted 
au(0: ^u(i)j etc.; this may be accomplished algebraically 
from the parameters of x(t), see Appendix E of Lilly 
and Gascard (2006). The ratio of the two quantities 
V{t) and R{t) defines a frequency ^{t) = V{t)/R{t) 
which we call the geometric frequency. 

A Lagrangian trajectory can then be represented 
as the sum of a number M of different oscillatory 
displacement signals x^™^(t), each of the form (3), plus 
a residual: 



x(t) 



x{t) 

y{t) 



M 

E 



E{™>(t)-|-e(t). 



(5) 
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The residual signal e{t) includes the turbulent back- 
ground flow, as well as any non-oscillatory component 
such as a mean flow or the systematic self-propagation 
tendency of a vortex. The oscillatory signals may be of 
finite duration and may be overlapping in time, so that 
zero, one, or more than one such signal may be present 
at each moment. 

The problem is then to estimate the oscillations 
xi"^^(t) given an observed trajectory x(t), and from 
these estimates to characterize the M different ellipse pa- 
rameters a{™>(t), 6{™>(i), 6'{™>(t), and so forth. 
This is accomplished using a method called "multivariate 
wavelet ridge analysis" {Lilly and Olhede, 2009a, 2011), 

leading to estimates x^™^ (t) of the M modulated oscilla- 
tions in each trajectory; here the subscript "-0" indicates 
the choice of filter or wavelet il^{t) used in the analysis. 
A discussion of the basic idea and implementation of the 
method, together with details of the parameter settings 
used here, are provided in a technical appendix. 

From the method, we obtain estimates of ellipse 
properties at each moment. As an example, observe the 
good agreement in Fig. 1 between the ellipses, inferred 
non-locally from individual particles, and the Eulerian 
structures in the potential vorticity field. 

4 Results 

The results of this analysis are shown in Figs. 2 and 3. 
Subtracting the sum over all M estimated oscillations 
in each time series x(i) leads to an estimate e(i) of the 
non-oscillatory background flow. Fig. 2b. Comparison 
with the original time series, Fig. 2a, shows that the 
tightly looping motions associated with vortices appear 
to have been nearly completely removed. What remains 
behind is observed to consist of curving but disorganized 
motions, together with systematic motion associated 
with the jet and with eddy drift. Even cycloidal features 
in Fig. 2a, typically low-frequency motion of a particle 
on the far flank of an eddy, are largely removed. The 
ability to perform such a separation on this relative 
large set of drifters, with no intervention by the analyst, 
by itself constitutes a technical breakthrough. 

Instantaneous ellipses associated with the estimated 
oscillatory signals are shown in Fig. 2c, d. Coloring the 
ellipses according to their degree of eccentricity and 
sense of rotation in Fig. 2c reveals nearly exclusively 
cyclonic motions (red) on the poleward side of the jet 
and anticyclonic motions (blue) on the equatorward 
side. The separation of vortices by their polarity in this 
system is due to nonlinear beta drift acting on eddies 
generated within the jet core during the initial adjust- 
ment. Close inspection reveals some very small circles 
of opposite color inside some of the eddies; as discussed 
in the technical appendix, these represent weak, low- 
frequency signals, likely associated with the presence 
of exterior opposing vorticity anomalies. Another type 
of variability, nearly linear in polarization and oriented 
meridionally, is observed along the jet axis (green) . 

A time scale distinction between these two different 
types of motions can be seen in Fig. 2d, where the 
color coding represents the instantaneous oscillation 
period 2t: /zuit) as deduced from the geometric fre- 
quency w{t). The highly eccentric motions along the 



jet axis are seen to be an order of magnitude lower 
in frequency than the vortex motions. The former 
arise as particles in the jet are swept eastward through 
meanders caused by the deflection of the jet axis by 
low-frequency Rossby wave variability, as is readily 
apparent in the animation file vortexmovie .mov at 
ftp : //f tp .nwra. com/outgoing/lilly/vortex . 

A more detailed view of the properties of the estimated 
elliptical signals is found in the distribution plots on 
the radius/velocity plane of Fig. 3. A histogram on the 
geometric radius R{t) / geometric velocity V{t) plane is 
formed in Fig. 3a by binning all time points associated 
with each estimated modulated oscillation in all 100 
trajectories. Note that the slope V{t)/R{t) on this plane 
is the geometric frequency wit). This histogram clearly 
reveals a Rankine-type vortex profile — a solid-body core 
plus l/r decay — for cyclonic motions, V{t) > 0. As may 
be expected, the sohd-body histogram peak hes along the 
slope corresponding to the maximum vorticity anomaly 
of the initial unstable jet profile (gray line); because 
under the assumption of solid-body rotation we have 
C — 2V/R, this occurs at a frequency of Co/2 = 0.40/. 

This vortex profile pattern emerges primarily through 
the superposition of a number of different Lagrangian 
particles at different locations in different eddies. On 
the anticyclonic side the pattern is less evident, since by 
chance, more particles have ended up in cyclonic eddies 
compared to anticyclonic eddies in this simulation; 
this asymmetry is a reminder of the well-known slow 
convergence of Lagrangian statistics owing to long-term 
particle trapping (e.g. Pasquero et al, 2002). 

The low-frequency motions are evident as the symmet- 
ric maxima in Fig. 3a around the horizontal line V = 0. 
The median eccentricity e in each bin. Fig. 3b, shows that 
the motions in the vortex profile are nearly circular in 
polarization, while the low-frequency motions are nearly 
linear, confirming that these two regions on the R/V 
plane correspond to the two types of features seen in 
Fig. 2c, d. An important point is that the nearly circular 
vortex motions and the highly eccentric low-frequency 
motions occupy almost non-overlapping regions of the 
radius/velocity plane, apart from very large-scale (say 
~50 km radius) and low-frequency motions, which 
may either be generated by jet meander or by circular 
oscillatory motions on the far flank of a vortex. 

5 Conclusions 

This paper has applied a new analysis method — 
multivariate wavelet ridge analysis — to the identification 
of oscillatory motions in Lagrangian trajectories from 
a numerical simulation of an unstable jet. It is shown 
that vortex motions can be effectively extracted from 
the trajectories and described locally in terms of their 
time-varying frequency content and ellipse geometry. 
Meridional meandering of the jet axis constitutes 
another strong oscillatory signal in this model, but 
these motions are clearly distinguished from the vortex 
motions on account of their much lower frequency, 
nearly linear shape polarization, and different location 
in radius/ velocity space. The method is therefore able 
to unravel the superposition of different processes in 
individual trajectories, making possible the investigation 
of these separate processes in isolation from one an- 
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Figure 2: Panel (a) shows all Lagrangian trajectories from the unstable barotropic jet simulation, dispersing from 
their initial location at y = 0, with each trajectory in a different color. In (b), the residual curves from the wavelet 
ridge analysis are presented. Snapshots of ellipses corresponding to the ellipse properties estimated from the 
wavelet ridge analysis are shown in (c) and (d). In both panels, highly eccentric ellipses with e > 0.95 — typically 
very low- frequency signals — are plotted with a time step of 1/2 the estimated period, while other ellipses are 
plotted with the time step of every 3 estimated periods. The highly eccentric ellipses are plotted in (c) in green, 
while the remaining anticyclonically rotating ellipses and cyclonically rotating ellipses are plotted in blue and red, 
respectively. Finally in (d), the ellipses are colored according to logj^p of the estimated geometric period 2TT/m{t), 
measured in days, as indicated in the color bar. 



other. These results serve as validation supporting the 
application of this method to large-scale observational 
studies. Not addressed here, but left to the future, are 
a detailed treatment of stochastic errors, the impact of 
measurement noise, comparison between the Eulerian 
and Lagrangian perspectives, and a consideration of 
method performance for other types of motions such as 
inertial oscillations or baroclinic instability waves. 
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supported by award #0751697 from the Physical 



Oceanography program of the United States National 
Science Foundation. The work of S. C. Olhede was sup- 
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Kingdom. The animation file was generated using 
Image2Movie by Jeffrey J. Early, available at http: 
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A Technical Appendix 

The extraction of individual oscillatory signals from 
the modeled Lagrangian trajectories is accomplished 
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Distiibutions of Estimated Oscillations in Unstable Jet Trajectories 




25 50 75 100 125 25 50 75 100 125 

Radius (km) Radius (km) 



Figure 3: Statistics of the estimated modulated elliptical signals shown on the geometric radius/velocity or R/V 
plane. Positive and negative velocities correspond to cyclonic and anticyclonic motions respectively. In (a) we 
see the histogram on the R/V plane occupied by all estimated oscillatory signals from all 100 trajectories, with 
a logarithmic color axis. The median eccentricity s in each bin is shown in (b). In both panels, the sloping gray 
lines correspond to ±Co/2, the geometric frequency value corresponding (for an eddy in solid-body rotation) to 
the maximum initial jet relative vorticity anomaly. The black lines are plus or minus the maximum and minimum 
cutoff frequencies in the analysis. 



using a method called "multivariate wavelet ridge 
analysis" {Lilly and Olhede, 2009a, 2011). All numer- 
ical code associated with this paper is distributed as 
part of a Matlab® toolbox called Jlab, available at 
http://www.jmlilly.net. The file makef igs_vortex 
provides the exact processing steps as well as scripts to 
generate all figures. This appendix describes the basic 
idea and presents an example. Equation numbers herein 
refer to the main text. 



An analytic wavelet is a time/frequency localized 
filter which has vanishing support on negative Fourier 
frequencies, and which is therefore complex-valued in 
the time domain, see e.g. Lilly and Olhede (2009b) and 
references therein for details. The wavelet transform 
of a real-valued signal vector x(t) with respect to the 



wavelet V'(i) is defined as 

/CO / ^ ^ \ 

and can be seen, on account of the 1/s normalization, 
as a set of bandpass operations indexed by the "scale" s 
which controls the dilation or contraction of the wavelet 
in time. The frequency-domain wavelet 4'(w) obtains 
a maximum magnitude, set to "^(uj) = 2, at some 
frequency uj^. This frequency is used to convert scale 
into a period via 2i:s/uj^. 

The particular choice of analyzing wavelet ^{t) is 
important in order to minimize errors in the subsequent 
analysis [Lilly and Olhede, 2009b, 2011). We use the 
"Airy wavelet" of Lilly and Olhede (2009b), which can 
be seen as a superior alternative to the popular but 
problematic Morlet wavelet, as discussed therein. The 
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Airy wavelet is controlled by a parameter P^i, , with -P,/, /tt 
giving the number of oscillations spanning the central 
window of the wavelet. We choose P^/n = ^/Q/tt ~ 0.78 
in order to obtain a high degree of time concentration 
at the expense of frequency resolution. 

Multivariate wavelet ridge analysis estimates modu- 
lated oscillations in a multivariate, or vector- valued, time 
series by first identifying maxima of the transform mag- 
nitude. A brief introduction to this method may be 
found in Lilly and Olhede (2009a), with further details 
and bias estimates provided by Lilly and Olhede (2011). 
A ridge point of Wx.^ {t, s) is defined to be a point on the 
{t,s) or "time/scale" plane satisfying 



transform evaluated along the ridge may be written as 



ds 



||wx,v,(i,s)|| =0, 



|wx,v-(i,s)|| < 



and thus ridge points are locations where the norm of the 
wavelet transform vector achieves a local maximum with 
respect to variations in scale. Adjacent ridge points are 
then connected to each other to yield a single-valued, 
continuous scale curve as a function of time called a 
ridge ?(t) . 

A number, say M, different ridges 's^'^^it) may be 
present in the same time series, and may overlap in time. 
The real part of the wavelet transform along the mth 
ridge may be taken as an estimate of the mth oscillatory 
signal x^"^^{t) in the composite model (5). That is 



j{m} 



(t) ^ 3? { wx.v, (t, s^™> it)) } = x{™> it) + a4"> (t) 



where the last term on the right-hand side is an error 
term. This may be understood as resulting from the 
fact that the real part of the wavelet transform at scale 
s = w^/w is a bandpass at frequency oj; thus the wavelet 
ridge estimate combines different pass bands at different 
times, as is appropriate for the time- varying nature of 
the signal. The error term depends upon several factors: 
(i) the local ratio of the oscillatory signal strength to 
the background variability strength, (ii) the magnitude 
of changes in oscillation properties over the wavelet 
timescale, see the treatment in Lilly and Olhede (2011), 
and (iii) the distance between contemporaneous oscilla- 
tory signals in frequency compared with the wavelet fre- 
quency profile. A detailed treatment of errors is outside 
the scope of this paper; instead, we will point to Fig. 2b 
in the text and the example below as subjective evidence 
that oscillatory signals have been recovered satisfactorily. 

It remains to form estimates of the ellipse parameters 
in the modulated elliptical signal model (3) associated 
with each of the M ridges. Observe that, as there are two 
parameters on the left-hand side of (3) but four param- 
eters on the right-hand side, the ellipse parameters are 
underdetermined for a given oscillatory signal x(t). Since 
the analytic wavelet transform Wx,^(t, s) is a complex- 
valued 2-vector, it consists of four quantities at each 
time. It has been shown by Lilly and Gascard (2006) and 
Lilly and Olhede (2010) that the complex-valued wavelet 



W: 



.,^(t,S^™>(i)) = 

[cos^{">(t) 
.sin^{™>(t) 



-sin^{'">(t)' 
cos^{">(t) . 



a{™>(t) 
-ife{'">(t) 



with the four quantities on the right-hand side being the 
sought-after estimates of the four ellipse parameters at 
each moment for ciach ridge. The "hats" mark estimated 
ellipse parameters that are defined implicitly through 
this equation. 

In practical implementation of the wavelet ridge 
analysis, two numerical thresholds must be introduced. 
Firstly a cutoff on the amplitude of ridge points is 
applied in order to reject very weak oscillations, which 
we choose here as position oscillations of less than 400 
meters (1/100 of the jet half- width Y) in magnitude. 
Secondly we define the ridge duration 



dt 



where "arg" denotes the complex phase, and the 
integral is taken over the time interval T^™^ of the mth 
ridge. The ridge duration gives the number of complete 
oscillations executed by the estimated signal; we reject 
all ridges that execute fewer than D^™} = 2P^,/tt w 1.6 
complete cycles, as we find ridges that are too short com- 
pared with the wavelet duration are generally spurious. 

An example of wavelet ridge analysis is presented in 
Fig. Al for one of the trajectories from the numerical 
simulation. The position trajectory x(i) shown in 
Fig. Ala is transformed using the parameter settings 
discussed above. 104 different scale levels are used which 
vary logarithmically from a minimum period 2t:s/uj^ of 
0.6 days to a maximum of 103 days. The norm of the 
wavelet transform vector is shown in Fig. Alb which, 
like position record itself, has imits of kilometers. In 
the wavelet transform we see low-frequency variability 
with periods of ~10 days, characteristic of geostrophic 
turbulence, and higher- frequency variability with peri- 
ods of ^1 day reflecting particle motion aroTind vortex 
cores. The large spectral gap between them supports 
our choice of a time- localized wavelet. 

M = 4 ridges are detected and are marked in the 
figure as curves indicating maxima of the wavelet trans- 
form with respect to scale. The real part of the wavelet 
transform along the mth ridge estimates the mth os- 
cillatory signal component in the composite model (5), 
which the analysis is able to isolate from the surrounding 
variability. These four estimated signals are then plotted 
in Fig. Ala along with the residuals. It is observed that 
the oscillatory features in the original trajectory are 
removed extremely well by the analysis method. 

In the central part of the record, two ridges are 
simultaneously present, with the lower-frequency curve 
representing an interaction with a distant vortex. A 
very weak amplitude ridge, just greater than our cutoff 
amplitude, is detected near the end of the time series; 
this is almost certainly spurious and illustrates the 
reason why we impose an amplitude cutoff; weak "false 
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Example of Wavelet Ridge Analysis 
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Figure 4: An illustration of wavelet ridge analysis. In (a), the east-west and north-south displacement signals 
associated with a Lagrangian "particle" in the model simulation are shown with the blue and red heavy lines, 
respectively. Thin blue and red lines show aggregate east-west and north-south oscillatory displacement signals 



obtained by summing over all estimated modulated oscillations (t) at each time. The green curves are the 

east-west and north-south residuals x(i) — X]m=i ^^'"'^(^)- The wavelet transform norm ||wx,^(i, s)|| is shown in 
(b), with scale converted into period on the y-axis. Four ridge curves are identified and drawn. The real part of 
the wavelet transform evaluated along these four curves give the estimated oscillations used in panel (a). Vertical 
lines in both panels mark the approximate times of vortex merger events as seen in the animation. 



positives" such as this signal do not impact our analysis 
because their amplitude is so small, and because they 
appear to occur rarely. Note that the low-frequency 
background, while slowly meandering in nature, does 
not generate ridges because it does not typically result 
in the appearance of multiple orbits through the same 
oscillatory structure, and is therefore not characterized 
as a modulated oscillation by this analysis. 

There are two transitions in this time series: one 
around 80 days, when the main ridge lowers its fre- 
quency and a second contemporaneous ridge appears, 
and a second transition around 125 days, when these 
two ridges disappear. It can be seen from watching the 
animation that these two transitions correspond to two 
vortex mergers. In that animation, ellipses associated 
with this trajectory are drawn in black, and the position 



of this drifter is indicated with a magenta dot. During 
the first merger, the vortex grows in size and its 
frequency decreases, but there is a remaining distant 
vortex which causes lower-frequency oscillatory motion 
that is manifested as the lower of the two ridges. During 
the second merger, the particle is ejected in a filament 
and consequently the oscillatory motions come to an 
end. This indicates that transitions in the ridges can 
have physical meaning and can capture vortex evolution. 
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